JOURNAL OF GEOPHYSICAL RESEARCH. VOL. 102, NO. D12, PAGES 13,555-13,569, JUNE 27, 1997 

NASA-CR-205091 

Tracer lamination in the stratosphere: A global 

climatology 

Christof Appenzeller 1 and James R. Holton 

Department of Atmospheric Sciences, University of Washington, Seattle 


Abstract. Vertical soundings of stratospheric ozone often exhibit laminated tracer structures 
characterized by strong vertical tracer gradients. The change in time of these gradients is used 
to define a tracer lamination rate. It is shown that this quantity can be calculated by the cross 
product of the horizontal temperature and horizontal tracer gradients. A climatology based on 
UARS satellite-borne ozone data and on ozone-like pseudotracer data is presented. Three strato- 
spheric regions with high lamination rates were found: the part of the stratospheric overworld 
which is influenced by the polar vortex, the part of the lowermost stratosphere which is 
influenced by the tropopause and a third region in the subtropical lower stratosphere mainly 
characterized with strong vertical shear. High lamination rates in the stratospheric overworld 
were absent during summer, whereas in the lowermost stratosphere high lamination rates were 
found year-round. This is consistent with the occurrence and seasonal variation of the hori- 
zontal tracer gradient and vertical shear necessary for tilting the tracer surfaces. During winter, 
high lamination rates associated with the stratospheric polar vortex are present down to - 100 
hPa. Several features of the derived climatology are roughly consistent with earlier balloon- 
borne studies. The patterns in the southern and northern hemisphere are comparable, but 
details differ as anticipated from a less disturbed and more symmetric southern polar vortex. 


1. Introduction 

The instantaneous distribution of a trace constituent in the 
atmosphere depends not only on its chemical sources and 
sinks but also on its redistribution due to transport induced by 
various dynamical processes. Large-scale advection such as 
vortex erosion and wave-breaking [e.g., McIntyre and Palmer , 
1984; Juckes and McIntyre , 1987; Hess and Holton , 1985; 
Hess , 1990; Bowman, 1993; Bowman and Mangus , 1993; 
Waugh et ai, 1994] together with differential vertical shear 
[e.g., Murphy et ai , 1989; Reid and Vaughan , 1991; Orsolini , 
1995] provide the potential for folding and tilting of tracer 
surfaces that lead to filament-like tracer structures in the 
stratosphere. Concerns about ozone depletion within the 
polar vortex, and subsequent mixing of ozone depleted air 
into lower latitudes, have triggered much interest in these 
processes [World Meteorological Organization ( WMO ), 1986] 
since the induced horizontal and vertical scale reduction 
implies the mixing of air masses from substantially different 
origins, both in a dynamical and chemical context [e.g., 
Murphy et ai, 1989; Pierce and Fairlie , 1993; Waugh et ai , 
1994; Plumb et ai , 1994; Chen, 1994; Haynes and Anglade, 
1997]. 

In a vertical sounding, tracer filaments can appear as layers 
of enhanced or depleted concentration referred to as laminated 
tracer structure. Ozone laminae of varying depth and magnitude 
have been observed since balloon soundings were available 
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[e.g., Dobson , 1973; Gardiner , 1988; Reid and Vaughan , 
1991] and have also been found in lidar observations [Reid 
et ai , 1993; Newman et ai, 1996; Orsolini et al. , 1997] and 
airborne in situ measurements [Murphy et ai , 1989]. 

The distribution of laminae is highly variable both in space 
and time. Due to the lack of global stratospheric ozone 
measurements, only a limited number of climatological studies 
have been performed. Amplitudes of the local variations of 
the ozone mixing ratio typically show a peak in the mid 
stratosphere around 10 hPa [I/.S. Standard Atmosphere , 1976; 
also Andrews et ai, 1987] and a seasonal variation with a 
maximum in midwinter and a minimum in midsummer [Roeth 
and Ehhalt, 1987]. A particular class of laminae characterized 
by an anomaly in ozone partial pressure [ Dobson , 1973; Reid 
and Vaughan , 1991; Reid et ai , 1993] has been observed to 
occur most commonly in winter and spring at high latitude and 
with a peak altitude of occurrence around 14 to 16 km. In the 
summer season such strong ozone laminae were absent. The 
preferred occurrence near the winter polar vortex suggests that 
the likely cause of ozone laminae is differential advection of 
air masses of different origin which are interleaved in the 
vertical. 

The above mentioned climatological studies are based on 
routinely performed balloon-borne ozone soundings using 
either a long time series for a single station or a pseudo- 
meridional cross section based on 13 sounding stations in 
the northern hemisphere and 5 in the southern hemisphere 
[Reid and Vaughan, 1991]. The recently available satellite- 
borne ozone measurements provide a much higher global 
coverage; however, these data have not yet been widely used to 
address the issue of tracer lamination. Satellite data typically 
have poor vertical resolution compared to balloon soundings. 
Fortunately, recently developed contour advection techniques 
provide a powerful tool to dramatically increase the effective 
vertical and horizontal data resolution [e.g., Waugh and 
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Figure 1. Schematic diagram illustrating the flow situation leading to high lamination rates. 


Plumb, 1994; Waugh et al ., 1994; Plumb et ai , 1994; Chen 
et al., 1994; Pierce et al., 1994; also Appenzeller et al . , 
1996a]. Satellite ozone data have been used to trace back air 
masses in which laminae are embedded [Orsolini et al. , 1995; 
Orsolini, 1995; Orsolini et al ., 1997; Newman and Schoeberl , 
1995; Newman et al., 1996] and have confirmed the 
involvement of air from the polar vortex. The successful 
contour advection simulations based on coarse wind analysis 
data suggest that comparatively large-scale flow features are 
responsible for the folding and tilting of tracer surfaces that 
finally lead to laminated structures. 

A characteristic feature of any lamina or filament-like 
structure found in a vertical tracer sounding is a strong anomaly 
in the vertical tracer gradient of either positive or negative 
sign. The development of a laminated tracer structure can 
therefore be characterized by a temporal change of the vertical 
tracer gradient. This quantity will be referred to as the tracer 
lamination rate and will be explored in this study. Two types 
of tracer data will be used: coarse resolution satellite ozone 
measurements and a more highly resolved pseudotracer data set 
based on the atmosphere’s potential vorticity (PV) structure. 
The results will be presented in the form of probability charts 
for the occurrence of high tracer lamination rates. This 
quantity can be interpreted as a proxy climatology of the 
production of laminated tracer (ozone) structures. 


2. Theory of Production of Vertical Tracer 
Gradients 

In Figure 1 the time development of a tracer lamina is 
schematically illustrated. The hypothetical tracer is assumed 
to have a horizontal distribution with a strong but localized 
tracer gradient near x — 0 and two different constant tracer 
concentrations outside this region. For such a distribution the 
surfaces (s) with constant tracer mixing ratio (%) are com- 
pletely vertically oriented (Figure lb) and a balloon sounding 
at any point in x would show a constant tracer concentration 
(Figure Id). Now assume that a vertical shear that is varying 
with height is acting on this tracer distribution (Figure la). 
After a certain time t the tracer surfaces will be tilted away 
from their initial positions, thereby creating vertical tracer 
gradients of positive and negative signs (Figure 1c). In a 
balloon sounding near x = 0, the tilted tracer surfaces now 
appear as a tracer filament (Figure le). The strongest filament 
occurs in conjunction with the strongest horizontal tracer 
gradients. In essence, differential vertical shear that is not 
aligned with the tracer surfaces converts horizontal tracer 
gradients into vertical ones. 

To quantify this mechanism, consider the conservation 
equation of a long-lived tracer with a volume mixing ratio 
(X ) : 



APPENZELLER AND HOLTON: TRACER LAMINATION CLIMATOLOGY 


13,557 


D 
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Here a log pressure coordinate system as in Andrews el al 
[1987] was used and 

D d d d d 
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represents the material derivative. S is the net source term 
and is assumed to be zero. Taking the vertical derivative of 
equation (1), the material rate of change of the vertical tracer 
gradient can be expressed as 

_£_f ^ | (2) 

Dt\dz J dz dx dz dy dz dz 


Using the geostrophic assumption the vertical shear is directly 
linked with the horizontal temperature gradient, i.e., 

fdu aA = JL?L) (3) 

\dz' dz) \ Hf dy' Hf dx) 

H represents the scale height,/ the Coriolis parameter and R 
the gas constant. Equation (2) can therefore be rewritten as 
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Since large-scale motions are quasi nondivergent, scaling 
suggests that the last term in equation (4) can be neglected 
compared to the term in the brackets and equation (4) simpli- 
fies to 


— k • [vr X V*] (5) 

Dt \ dz y Hf 

The magnitude of the cross product of the horizontal gradients 
(VT, V*) is therefore a measure of the production of vertical 
tracer gradients following the motion. The absolute value of 
this production rate will be referred to as the tracer lamination 
rate (£): 
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The three-dimensional distribution of £ can be calculated from 
isobaric temperature and tracer distributions. This is particu- 
larly useful for satellite-measured tracer data since these data 
have a comparatively low vertical resolution. For the analyses 
shown below, equation (6) was evaluated using spherical 
coordinates. 

In the winter stratosphere near the edge of the polar vortex, 
typical horizontal temperature gradients are V7 1 «20 x 1CL 6 
K m -1 and typical ozone gradients are V# ~1.5 x 10" 6 ppmv 
m _1 . For a hypothetical flow situation where the vertical 
shear is acting perpendicular to the tracer surfaces these values 


would result in lamination rates up to £ =10x 10 9 ppmv 
m _1 s _1 (with R/Hf « 400 J Kr 1 kg" 1 m _1 s). Or equivalently, 
to produce a hypothetical lamina with a relative amplitude of 
- 0.5 ppmv and a vertical depth of ~ 2 km about half a day 
would be needed. 

3. Data 

3.1. MLS Ozone and Temperature Data 

The microwave limb sounder (MLS) from the Upper 
Atmosphere Research Satellite (UARS) provides stratospheric 
ozone and temperature measurements that can be used to 
construct daily global isobaric ozone charts (for more details 
about the MLS data see, e.g., the special issue dealing with 
UARS result in Geophysical Research Letters , 20 , (12), 1993). 
In this study MLS Level 3b data were used and the 13 Fourier 
components were gridded onto a coarse longitude latitude grid 
(with a lati-tudinal coverage between 80°S and 80°N and a 
horizontal resolution of 24° longitude and 4° latitude). In the 
vertical (i.e., between 100 and 1 hPa) the data have been 
interpolated onto UARS standard pressure levels. For the mid 
to upper stratosphere the data are of high accuracy [Froidevaux 
et al . , 1996] but the data quality decreases for pressure greater 
than ~50 hPa. Therefore in the lowest stratosphere the ozone 
data should not be used for any physical interpretation. The 3 
year averages shown below were constructed using all 
available data in the period from November 1991 to October 
1994. 

3.2. Pseudotracer Potential Vorticity 

The second data set used is a pseudotracer based on the atmo- 
sphere’s potential vorticity (PV) structure. Since on synoptic 
timescales both PV and ozone are roughly conserved quantities 
following an air parcel and both tracers have a comparatively 



Ozone (ppmv), MPV (pvu, scaled by 2/3) 


Figure 2. Vertical distribution of ozone (solid curve, given 
in parts per million by volume) and modified potential vorti- 
city MPV (dashed curve, given in potential vorticity units 
scaled by 2/3) for midlatitude winter based on U.S. standard 
atmosphere data. Dashed-dotted curve is MPV based on UKMO 
January climatology for 1993 at latitude 45°N. 
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strong meridional gradient, it can be anticipated that there 
exists a strong similarity in the pattern of their horizontal 
distribution. Strong horizontal gradients in both tracers can 
in particular be anticipated near the polar vortex edge and the 
tropopause. Unlike ozone, Ertel’s PV as conventionally defined 
increases exponentially with height. Scaling of the Ertel PV 
as suggested by Lait [1994], 



allows for a vertical PV profile more comparable to that of 
ozone (Figure 2) without losing the fundamental conserva- 
tion properties [Ertel, 1942]. The modified potential vorti- 
city (MPV) was calculated from daily (00 GMT) global 
stratospheric analyses prepared by the United Kingdom 
Meteorological Office (UKMO) [Swinbank and O'Neill , 1994]. 
Horizontal winds, geopotential heights, and temperature data 
were used with a horizontal resolution of 2.5° latitude by 3.75° 
longitude. In the vertical, UARS standard pressure levels were 
used corresponding to a resolution of roughly 2.5 km ranging 
from - 1000 hPa to 1 hPa. The pseudotracer is displayed 
in standard PV units with 1 pvu = 10 -6 Km 2 kg -1 s _I . 
Pscudotracer lamination rates were not calculated at the two 
lowest pressure levels (i.e., 1000 and 681 hPa) as well as near 
the equator (i.e., between -10 and 10 latitude) and at the poles. 
It is important to keep in mind that these data represent a 
blend of sparsely available operational stratospheric obser- 
vations with the model’s own 6 hour forecast. Some of 
the small-scale PV features might therefore be artificially 
generated by the analysis procedure. 

4. Results 

4.1. Maps 

Plate la shows an example of the northern hemisphere 
stratospheric polar vortex at 10 hPa as seen by the MLS 
measurements. Ozone (black curves) and temperature (color) 
are shown. The temperature distribution was dominated by a 
wave number 2 pattern with cold temperatures over the Atlantic 
and warm ones over Siberia providing several regions with 
strong temperature gradients. The ozone pattern revealed that 
the polar vortex was shifted toward northern Europe and was 
characterized by relative low ozone concentrations compared 
to more southern latitudes. Relatively strong ozone gradients 
were found at the edge of the polar vortex and also at 
midlatitudes stretched across the Atlantic to the Pacific. As 
anticipated, the overall horizontal pattern of the pseudotracer 
MPV truncated to the low MLS resolution (Plate 2a) closely 
resembled the observed ozone distribution (Plate la). Note the 
opposite sign in meridional gradient. 

In Plates lb and 2b the corresponding tracer lamination 
rates (equation (6)) are shown. Three regions with high lami- 
nation rates were found for both tracers, all of which are linked 
to the edge of the vortex: one over Europe, one over Siberia 
and a third over northern Canada. In these regions the ozone 
(pseudotracer MPV) gradients were relatively strong and not 
parallel to the temperature gradients as anticipated. Over the 
midlatitude Atlantic and the Pacific, the lamination rates were 
substantially lower despite relatively strong ozone (MPV) 
gradients. 

The same calculations are shown in Plate 3, but using the 
pseudotracer with full UKMO data resolution. The higher- 


resolution data set leads to a much finer-scaled and more 
localized tracer lamination field. Spots with high values 
occurred again in the three regions mentioned above and were 
in conjunction with strong pseudotracer gradients. Other 
small-scale features like the one over eastern Siberia (160° E, 
60° N) lead also to enhanced production rates. These pseudo- 
tracer lamination rates were typically a factor 3/2 higher than 
for the low-resolution data and a factor 10 higher compared to 
low-resolution ozone. A second illustrative example in the 
lower stratosphere (-100 hPa) is given in Plate 4. A tongue of 
high pseudotracer values extends southwestward from of the 
polar vortex stretching across Eurasia into a region with high 
temperature. Enhanced absolute values of tracer lamination 
rates occurred on both sides of the streamer, indicating that the 
vertical tracer gradients would increase on the southern and 
decrease on the northern side. 

4.2. Vertical Profiles 

Two vertical MLS ozone soundings taken within the region 
of high lamination rates over Europe (Plate lb) are shown in 
Figure 3. The vertical profiles above 10 hPa showed stronger 
vertical gradients compared to a 3 year monthly mean taken at 
the same grid points. This is consistent with the enhanced 
lamination rates shown above. Due to the poor horizontal 
and vertical resolution of the satellite data ( i.e., roughly one 
measurement every 5 km in the vertical) the soundings can 
only resolve comparatively large-scale tracer structures and 
any smaller-scale laminae-like features are not resolved. 

4.3. Proxy Climatology 

In order to construct a global and seasonal climatology the 
tracer lamination rates were calculated for each day that tracer 
data were available (between November 1991 and October 
1994). Subsequently the values were grouped into two main 
categories (see Table 1), group 1 containing no to weak 
lamination rates, group 2 moderate to strong lamination rates. 
A third group (group 3) was introduced to have some control 
over extremely high values. The combined group 2 and 3 will 
simply be referred to as high lamination rates. The results are 
given in the form of averaged probability charts (ranging from 

0 to 1) to observe a particular group of lamination rates. Note 
that the longitude-latitude grid used is not equidistant in kilo- 
meters. The threshold values were chosen subjectively based 
on inspection of daily charts (compare with Plates 1 to 3) and 
were kept as simple as possible. Moderately changing the 
threshold value did not change the overall conclusion given 
below. 

Meridional cross sections of monthly averaged charts of the 
probability to observe high (combined group 2 and 3) lamina- 
tion rates are shown for low-resolution MLS ozone (Figure 4) 
in the altitude range 100 to 1 hPa and for the high-resolution 
pseudotracer MPV (Figure 5) in the altitude range 1000 to 

1 hPa. The data represent a 3 year average for the time periods 
January-February (Figures 4a and 5a) and July-August (Figures 
4b and 5b). 

In the stratospheric overworld a clear difference between 
winter and summer season was found. In the winter season at 
midlatitudes and high latitudes high lamination rates occurred 
in both the MLS ozone ( -20 to 30% of all the cases) and in 
pseudotracer MPV (-30 to 40% of all the cases). In the summer 
season, even moderate production rates were absent. The 
pattern further suggests a slight asymmetry between the two 
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Plate 1 . (a) Northern hemisphere ozone and temperature distribution as seen by MLS at 10 hPa on 

March 18, 1992. Ozone volume mixing ratio is given by black curves (black dotted curves below 8 ppmv, 
interval 0.5 ppmv) and temperature is given in color (contour interval 2.5 K). The area shown is northward of 
20°N. (b) Same as Plate la, but calculated ozone lamination rates (colored, values given in units of 10" 9 ppmv 

m _1 s _1 , contour interval 0.4 x 10“ 9 ppmv m _1 s _1 ) and ozone volume mixing ratio overlaid (contours as in 
Plate la). 
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Plate 2. (a) Same as Plate 1 but for pseudotracer MPV based on truncated UKMO potential vorticity 

structures. MPV is given by black curves (dotted curves below 10 pvu, interval 2 pvu) and temperature in color 
(contour interval 2.5 K). (b) Same as Plate 2a, but pseudotracer lamination rates based on truncated data 
(colored, values given in units of 10“ 9 pvu m _1 s _1 , contour interval 2.5 x 10 -9 ppmv m _1 s _l ) and pseudotracer 
MPV overlaid (contours as in Plate 2a). 
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Plate 3. Same as Plate 2b, but pseudotracer MPV lamination rates based on full data resolution (colored, 
contour interval 2 x 10 -9 pvu nHs -1 ) and pseudotracer MPV overlaid (contour interval 2 pvu). 


hemispheres. In the northern hemisphere the number of events 
increases with latitude, whereas in the southern hemisphere 
the region where high production rates occur seems to be more 
localized around a maximum at -60° south. 

The analyses based on the pseudotracer (Figure 5) suggest 
two other regions in the atmosphere where high lamination 
rates occur, the lowermost stratosphere and the subtropical 


lower stratosphere (in both these regions no reliable MLS 
ozone data were available). In the midlatitudes, high values 
typically occurred between 400 and 200 hPa, whereas in the 
subtropical region the maximum is located between -120 hPa 
and 60 hPa. In contrast to the variation in the stratospheric 
overworld the high lamination rates in the lowermost 
stratosphere were also present in the summer months. 



Plate 4. Lamination rates for pseudotracer MPV at 100 hPa over Eurasia for January 1, 1993. Lamination 
rates are colored (contour interval 2 x 1C T 9 pvu nrV 1 ), pseudotracer MPV overlaid in solid contours (contour 
interval 1 pvu), and temperature overlaid in dashed contours (contour interval 2.5 K). 
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Figure 3. Vertical ozone distribution as seen by MLS at grid 
points (12°E, 56°N) (left curves) and (36°E, 56°N) (right curves) 
for March 18, 1992 (solid curves) and for a monthly mean 
(3 year average) at the same grid points (dashed curves). 
Locations are marked by stars in Plate lb. Values given in 
parts per million by volume and curves shifted by 4 ppmv. 


The complete seasonal variation of the probability of 
observing high lamination rates is shown in latitude time 
sections in Figure 6 (for MLS ozone) and Figure 7 (for the 
pseudotracer MPV). For better display the seasonal cycle is 


Table 1 . Threshold Values Used for Grouping the Tracer 
Lamination Rates (£) 


Group 

Description 

No. 

£ MLS Ozone £ Pseudotracer 
MPV 

Weak 

lamination 

no to weak 

1 

£ < 1 x 10“ 9 £ < 10 x 10~ 9 

Strong 

lamination 

moderate to 
strong 

2 

1 x 10- 9 10 x 10" 9 

< £ < 3 x 10 -9 < £ <30 x 10~ 9 


extremely 

strong 

3 

£>3xl0- 9 £>30xl0- 9 


Units of £ are parts per million by volume for ozone and 
potential vorticity units per meter per second for pseudotracer 
MPV. 


shown twice. The data have been vertically averaged in three 
subregions of the stratosphere. 

In the stratospheric overworld (Figures 6 and 7a), high 
lamination rates for both tracers occurred during the winter 
half-year and mainly poleward of 30°. For the southern 
hemisphere the time period was longer than for the northern 
hemisphere, consistent with the longer southern hemisphere 
winter. The period with high lamination rates was slightly 
longer in MLS ozone than for the pseudotracer data. In the 
northern hemisphere the highest probability was again found 


a) 



0.5 

0.4 

0.3 

0.2 

0.1 
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latitude (degree) 

Figure 4. Latitude-height section of the probability of finding high ozone lamination rates for the period 
(a) January-February and (b) July-August. Data are based on MLS measurements and averaged over 3 years 
(1992 to 1994). Vertical axis ranges from 100 to 1 hPa. Data at 100 hPa and near the equator (+/-10 0 latitude) 
are omitted. Contour interval is 0.025. 
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Figure 5. Same as Figure 4 but for the probability of finding high pseudotracer lamination. Data are based 
on UKMO analyses. Vertical axis ranges from 1000 to 1 hPa. Data below 650 hPa and near equator are omitted. 
Contour interval is 0.025. 


at high latitudes, whereas in the southern hemisphere the 
maximum occurs again around 60°S (midwinter) and moved 
poleward during the spring season. Such a distribution is 
consistent with a less disturbed and more symmetric southern 
hemisphere polar vortex. 

In the lowermost stratosphere (Figure 7b), high pseudo- 
tracer lamination rates were present throughout the year, but 
with substantially less probability and less horizontal extent 
during the summer season. The seasonal variation of the lati- 
tudinal coverage is consistent with the seasonal variation of 
the tropopause position [Appenzeller et al , 1996b]. In both 
hemispheres during the winter season a maximum occurred 
around +/-30 0 to +/-40 0 latitude near the subtropical tropo- 
pause break. It is interesting to note that Dobson [1973] 
suggested this to be the region where ozone poor air is 
entering the stratosphere. 

Similar features can be observed at 100 hPa (Figure 7c), but 
at higher latitudes the lamination rates showed a more pro- 


nounced seasonal cycle. In the northern hemisphere, high 
tracer lamination rates occurred only from January to mid-May 
and were completely absent from late spring to fall. A 
comparable seasonal cycle was found in earlier balloon studies 
[Reid and Vaughan , 1991] that examined the frequency of 
partial pres-sure laminae at the same altitude. These high 
lamination rates are associated with the lower edge of the 
stratospheric polar vortex. In the southern hemisphere the 
corresponding high probability region was observed below 
100 hPa, i.e., clearly in the lowermost stratosphere (Figure 
7b). At 100 hPa the location of the maximum moved again 
poleward during the transition from winter to spring. In the 
northern hemisphere only a very weak indication of such a 
pattern could be found. 

Finally, we note that extremely high lamination rates for 
both ozone and pseudotracer (characterized by the extra group 3 
in Table 1) were only found in the stratospheric overworld (not 
shown separately). They occur in the 3 winter months pole- 
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Figure 6. Seasonal variation of the probability of finding high ozone lamination rates in the stratospheric 
overworld. Data as in Figure 4, but a zonal and vertical (between 68 and 1 hPa) average. For display purposes 
the 3 year average is shown twice. Contour interval is 0.025. 


a) 



latitude (degree) 

Figure 7. Same as Figure 6, but seasonal variation of the probability of observing high pseudotracer 
lamination rates. Data as in Figure 6. (a) Probability vertically averaged in the stratospheric overworld 

(between 68 and 1 hPa). (b) Probability vertically averaged in the lowermost stratosphere (between 316 and 
146 hPa). (c) Probability at 100 hPa. Contour interval is 0.025. 
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latitude (degree) 
Figure 7. (continued) 


ward of ~60°N and in the southern hemisphere around 60°S. 
The maximum probability of observing a class 3 lamination 
rate was -10-15% in the northern and -5-10% in the southern 
hemisphere. At 100 hPa and in the lowermost stratosphere 
their contribution was negligible. 

5. Discussion 

Several features of the derived climatology are consistent 
with earlier studies based on balloon ozone soundings. These 
features include the preferred occurrence of filaments near the 


edge of the polar vortex [Reid and Vaughan, 1991; Reid et ai , 
1993, 1994], the preferred occurrence near the subtropical and 
midlatitude tropopause [ Dobson , 1973], an increase in tracer 
variability in the mid stratosphere [Roeth and Ehhalt, 1987; 
U.S. Standard Atmosphere , 1976], a strong seasonal cycle 
with weak variability in the summer season [Dobson, 1973; 
Reid and Vaughan , 1991; Roeth and Ehhalt , 1987], and a shift 
of the maximum of occurrence toward the pole during the tran- 
sition from winter to spring [Reid et ai , 1993]. Other features 
found by the balloon studies, such as the absence of laminae in 
the mid to upper stratosphere and in the subtropical lower 
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latitude (degree) 

Figure 8. Zonally averaged magnitude of horizontal gradients in (a) pseudotracer MPV and (b) UKMO 
temperature. Data are averaged over 3 years (1992 to 1994) for the period January-Fcbruary. Vertical axis 
ranges from 1000 to 1 hPa. (a) Values given in 10‘ 6 pvu m _1 and contour interval 1 x 10' 6 pvu m _1 . (b) Values 
given in 10* 6 K m ] and contour interval 1 x 10' 6 K m“*. 


stratosphere [ Dobson , 1973; Reid and Vaughan , 1991], as 
well as the occurrence of laminae at high latitudes well into 
early June, are not found in the derived climatology. Part of 
this difference is due to the different laminae concept applied 
(see below). 

Most of the distribution of high lamination rates can be 
understood in terms of the occurrence of both strong hori- 
zontal tracer gradients and the dynamical activity providing 
the vertical shear for tilting the tracer surfaces. Figure 8a 
shows a latitude-height section of the zonally averaged mag- 
nitude of horizontal pseudotracer gradients. As anticipated, 
regions with high averaged horizontal tracer gradients such as 
the tropopause region and the polar vortex region coincide 
roughly with regions characterized by high lamination rates 
(Figure 5). Regions with strong zonally averaged magnitude 
in vertical shear are shown in Figure 8b (following equation 
(3), the shear is represented by the horizontal temperature 


gradient). High shear values also exist in regions with low 
lamination rates. An example is the subtropical lower strato- 
sphere where strong vertical shears are present in both hemi- 
spheres but high pseudotracer lamination rates were only 
found in the northern (winter) hemisphere. 

The summer middle stratosphere is not only characterized by 
weak pseudotraccr gradients but also by weak vertical shear. 
The latter is consistent with the seasonal variation in wave 
activity in the stratospheric overworld that is strongly deter- 
mined by whether the stratospheric flow allows upward propa- 
gation of synoptic scale tropospheric wave sources Linear 
wave theory [Charney and Drazin , 1961; Andrews et al. y 
1987] indicates that upward propagation of Rossby waves is 
inhibited in easterly winds. Figure 9 shows the seasonal 
variation of the probability of observing a west wind (i.e., 
allowing upward propagation) in the stratospheric overworld. 
Regions with probability higher than 80% are shaded. The 
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Figure 9. Seasonal variation of the probability of observing 
a west wind in the stratospheric overworld (averaged between 
68 and 1 hPa). Regions with probabilities higher than 0.8 are 
solid shading, probability 0.6 is shown by solid lines, and 
probabilities less than 0.6 are dashed lines at 0.2 interval. 


are comparable or smaller [ Duetsch , 1978] than a typical 
threshold value used and therefore laminae defined on partial 
pressure anomalies are difficult to detect. 

The observed climatology in the lower and lowermost 
stratosphere relies solely on the pseudotracer MPV. MLS 
ozone data could not be used to verify the existence of lami- 
nated ozone structures below the 50 hPa surface. The pseudo- 
tracer has a somewhat better resolution than MLS ozone, but 
its distribution has only an approximate relation with that of 
ozone. For example, its maximum occurs at a lower altitude 
than the ozone maximum, and it decays faster in the upper 
stratosphere than does ozone (see Figure 2). However, accord- 
ing to equation (6) the tracer lamination rates are proportional 
to the absolute value of the horizontal tracer gradients and not 
to the value of the tracer concentration itself. For the tropo- 
pause region both tracers are characterized by comparatively 
strong gradients [e.g., WMO , 1986], whereas for the region in 
the subtropical lower stratosphere only a little information is 
available that supports the existence of strong horizontal 
ozone gradients [ Duetsch , 1978; Fahey et ai, 1996]. Finally, 
note that the occurrence of high lamination rates within the 
winter polar stratospheric vortex might be overestimated by 
small-scale structures artificially generated by the analysis 
scheme. Ozone data with better horizontal resolution cover- 
ing the entire stratosphere would substantially improve the 
quality of this analysis. 


pattern reflects several features found in the seasonal variation 
of strong tracer lamination rates. Upward propagation is also 
limited to the winter half-year and occurs only poleward of 
approximately 30° latitude. Again the winter period is longer 
in the southern hemisphere than in the northern, which is 
consistent with the less disturbed southern hemisphere’s polar 
vortex. 

These observations suggest that in the lower and lowermost 
stratosphere the occurrence of strong tracer lamination rates is 
mainly controlled by the existence of strong horizontal tracer 
gradients, since vertical shear is present throughout the year. 
In the mid to upper stratosphere, however, both the vertical 
shear and the strong horizontal gradients show a distinct 
seasonal cycle. 

The climatological patterns of the occurrence of laminated 
ozone structures are dependent on how exactly a tracer or 
ozone lamina is defined. In the current study the occurrence of 
a high lamination rate that measures a strong increase or 
decrease of vertical gradient in tracer mixing ratio is used as a 
proxy for the occurrence of a laminated tracer structure. No 
limitation is applied on its vertical scale. The lamination rate 
itself is inferred from the tilting of the large-scale horizontal 
tracer gradient by the large-scale vertical shear. The balloon 
studies [ Dobson , 1973; Reid and Vaughan , 1991; Reid et ai, 
1993] used a laminae definition based on a fixed ozone partial 
pressure {p 0 ) anomaly of at least 2 x 10" 3 Pa and a vertical 
extent of less then 2.5 km. Since ozone partial pressure 
decreases exponentially with height (i.e., p = v Q with 
the volume mixing ratio of ozone and p the atmospheric 
pressure), a laminae climatology based on a fixed partial 
pressure anomaly substantially underestimates the variations 
found in the stratospheric overworld. The same is true for 
laminae frequency found in the subtropical lower stratosphere, 
since in this region the absolute ozone partial pressure values 


6. Conclusion 

Vertical soundings of stratospheric ozone often exhibit 
laminated tracer structures characterized by strong vertical 
gradients in mixing ratio. These laminae indicate the possible 
mixing of air masses of dynamically and chemically substan- 
tially different origin. Examples are the potential mass trans- 
port across the relatively isolated stratospheric polar vortex 
[e.g., Murphy et ai , 1989; Waugh et ai, 1994; Chen et al . , 
1994] or across the extratropical tropopause [e.g., Holton 
et ai, 1995]. Both of these mixing processes substantially 
influence the atmosphere’s ozone distribution and chemistry. 
Dynamically induced variations in stratospheric tracer concen- 
tration are also important for data assimilation purposes. 
Since most of the laminae are at or below the current model 
resolution, tracer assimilation in regions and seasons with 
high tracer lamination may prove to be an extremely difficult 
task. 

The aim of this study was to provide a global and seasonal 
climatology of laminated ozone-like tracer structures in the 
stratosphere. Using the geostrophic assumption it was shown 
that production rates of vertical tracer gradients can be calcu- 
lated by the cross product of the isobaric temperature and tracer 
gradients. This method permits calculation of tracer lamina- 
tion rates based on satellite data with relatively low vertical 
resolution. The climatology presented was based on grouping 
lamination rates calculated from daily three-dimensional data 
sets rather then counting explicitly observed tracer laminae as 
in earlier studies. A high tracer lamination rate was taken as a 
proxy for the production of a laminated tracer structure. Two 
tracer data sets have been used, UARS satellite-borne low- 
resolution ozone measurements and higher-resolution ozone- 
like pseudotracer data based on modified potential vorticity. 
The ozone data were only available in the stratospheric over- 
world; the pseudotracer covered the entire stratosphere and was 
calculated from UKMO assimilated analyses. Comparisons of 
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low-resolution charts of the corresponding lamination rates 
showed a qualitative good agreement. A similar calculation 
with more highly resolved pseudotracer data showed more 
localized regions with higher production rates. Comparison 
with single point ozone soundings would require a much 
higher-resolution data set. 

The climatological analyses based on ozone and ozone-like 
pseudotracer data suggest that the stratosphere shows basi- 
cally three regions with frequent occurrence of high ozone 
lamination rates. These are the part of the stratospheric over- 
world which is influenced by the polar vortex, the part of the 
lowermost stratosphere which is influenced by the tropopause 
and a third region in the subtropical lower stratosphere. Both 
the polar vortex and the tropopause are characterized by 
relatively strong horizontal ozone gradients. High lamina- 
tion rates in the stratospheric overworld occur only during the 
winter half-year and are absent during summer, whereas high 
lamination rates in the lowermost stratosphere were found 
year-round. This is consistent with the occurrence and 
seasonal variation of the horizontal tracer gradient and verti- 
cal shear necessary for tilting the trace surfaces. The patterns 
in the southern and northern hemispheres are comparable, but 
details differ as anticipated from a less disturbed and more 
symmetric southern polar vortex. In the subtropical lower 
stratosphere a third region with high pseudotracer lamination 
rates was found. This region coincides mainly with strong 
vertical shear. Its relevance for ozone lamination is unclear. 
High lamination rates associated with the stratospheric polar 
vortex are present down to -100 hPa. In the northern hemi- 
sphere these high lamination rates occurred only from January 
to mid-May, which is roughly consistent with the seasonal 
cycle found in balloon-borne ozone partial pressure laminae. 
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